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We have implemented the new IR-improved Dokshitzer-Gribov-Lipatov-AltareUi-Parisi-Callan-Symanzik 
(DGLAP-CS) kernels recently developed by one of us in the HERWIG6.5 environment to generate a new 
MC, HERWIRI1.0(31), for hadron-hadron scattering at high energies. We present MC data that illustrate the 
comparison between the parton shower generated by the standard DGLAP-CS kernels and that generated by 
the new IR-improved DGLAP-CS kernels. We also show comparisons with FNAL data and we discuss possible 
implications for LHC phenomenology. 



1. Introduction 

In the LHC era of precision QCD, which entails 
predictions for QCD processes at the total preci- 
sion [1] tag ^ of 1% or better, we need resummed 
©(a^L"), 0{a,aL"-'), 0{a^L""),n = 0,1,2, n' = 
0, 1, 2, n" — 1,2 corrections, in the presence of parton 
showers, on an event-by-event basis, without double 
counting and with exact phase space. The roles of 
QED and EW effects [2, 3] are integral parts of the 
discussion with which we deal by the simultaneous 
resummation of QED and QCD large infrared(IR) ef- 
fects, QED (g) QCD resummation [4] in the presence 
of parton showers, to be realized on an event-by-event 
basis by MC methods; for, as shown in Refs. [3], no 
precision prediction for a hard LHC process at the 1% 
level can be complete without taking the large EW 
corrections into account. 

In what follows, we first review our approach to re- 
summation and its relationship to those in Refs. [5, 6]. 
Section 3 contains a summary of the attendant new 
IR-improved DGLAP-CS [7, 8] theory [9, 10]. Sec- 
tion 4 presents the implementation of the new IR- 
improved kernels in the framework of HERWIG6.5 [11] 
to arrive at the new, IR-improved parton shower MC 
HERWIRIl.O. We illustrate the effects of the IR- 
improvement first with the generic 2— >-2 processes at 
LHC energies and then with the specific single Z pro- 
duction process at LHC energies. We compare with 
recent data from FNAL to make direct contact with 
observation. Section 5 contains our summary remarks. 

For reference purposes, we call attention to the 
analyses in Refs. [12, 13], wherein the authors have 
argued that the current state-of-the-art theoretical 



precision tag on single Z production at the LHC is 
(4.1 ± 0.3)% (1.51 ± 0.75)%{QCD) ® 3.79{PDF) © 
0.38 ± 0.26{EW)% and that the analogous estimate 
for single W production is ~ 5.7%. One cannot em- 
phasize too much that these estimates show how much 
work is still needed to achieve the desired 1.0% total 
precision tag on the two respective processes, for ex- 
ample. 



2. QED^QCD Resummation 

We make use of the discussion in Refs. [4, 9, 10], 
wherein we have derived the following expression for 
the hard cross sections in the SM SU2L x [/i x SU!^ 
EW-QCD theory 



,~ _ SUMinXQCED) v^°o 1 r TT" 

"Oexp — e Z^n,m=0 n!m! J -l-'-Jl=l 



J2 = l J J2^^ 

Pn,m {kl, ■ ■ ■ , kn', kl, ■ ■ ■ , fcm) 0^ 0^ i (1) 



^By total precision of a theoretical prediction we mean the 
technical and physical precisions combined in quadrature or 
otherwise as appropriate. 



where the new YFS-style [14] residuals 

(3n,m{ki, . . . ,kn;k[, . . . ,k'^) have n hard gluons 
and m hard photons and we show the final state 
with two hard final partons with momenta p2, q2 
specified for a generic 2/ final state for definiteness. 
The infrared functions SUMir(QCED), Dqceb 
are defined in Refs. [4, 9, 10]. This is the exact, 
simultaneous resummation of QED and QCD large 
IR effects. 

One can see that our approach to QCD resumma- 
tion is fully consistent with that of Refs. [5, 6] as fol- 
lows. First, Ref. [15] has shown that the latter two 
approaches are equivalent. By using the color-spin 
density matrix realization of our residuals, we show in 
Refs. [9, 10] that our approach is consistent with that 
of Refs. [5] by exhibiting the transformation prescrip- 
tion from the resummation formula for the theory in 
Refs. [5] for the generic 2 — n parton process as given 
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in Ref. [16] to our theory as given for QCD by re- 
stricting (1) to its QCD component. In this way, we 
capture the respective full quantum mechanical color- 
spin correlations in the results in Ref. [16]. 



3. 1 R-lm proved DGLAP-CS Theory 

We show in Refs. [9, 10] that the result (1) allows us 
to improve in the IR regime ^ the kernels in DGLAP- 
CS [7, 8] theory as follows, using a standard notation: 



-l + z' 



{i-zy 



-/,(7,)5(l-^)], 



l,,l + (l-£f 



z 

l-z 



-/g(7g)<5(1-2)}, 

^G^'W = FYFs{jG)ei'''l{z'{l - zy^ + (1 - zfz'^'^}, 

(2) 

where the superscript "cxp" indicates that the ker- 
nel has been resummed as predicted by (1) when it 
is restricted to QCD alone - see Refs. [9, 10] for the 
corresponding details. These results have been imple- 
mented by MC methods as we exhibit in what follows. 

Let us first note that a number of ilhistrativc re- 
sults and implications of the new kernels have been 
presented in Refs. [9, 10, 19]. Here, we call attention 
to the new scheme [10] which we now have for preci- 
sion LHC theory: in an obvious notation, 

0- = ^ dxidx2Fi{xi)Fj{x2)d-{xiX2s) 
id 

= / dxidx2F'i{xi)F'j{x2)a-'{xiX2s), 



where the primed quantities are associated with (2) 
in the standard QCD factorization calculus. We 
have [4] an attendant showcr/ME matching scheme, 
wherein, for example, in combining (1) with HER- 
WIG [11], PYTHIA [27], MC@NLO [28] or new 
shower MC's [29], we may use either pr-matching or 
shower-subtracted residuals 

{0n,m{h,---,kn;k[,...,k'^)} to create a paradigm 



^This should be distinguished from the also important re- 
summation in parton density evolution for the "z 0" regime, 
where Regge asymptotics obtain - see for example Ref. [17, 18]. 
This improvement must also be taJsen into account for precision 
LHC predictions. 



without double counting that can be systematically 
improved order-by order in perturbation theory - see 
Refs. [4]. 

The stage is set for the full MC implementation of 

our approach. We turn next to the initial stage of this 
implementation - that of the kernels in (2). 



4. MC Realization of IR-lmproved 
DGLAP-CS Theory 

In this section we describe the implementation of 
the new IR-improved kernels in the HERWIG6.5 en- 
vironment, which results in a new MC, which we de- 
note by HERWIRIl.O, which stands for "high energy 
radiation with IR improvement""^. 

Specifically, our approach can be summarized as 
follows. We modify the kernels in the HERWIG6.5 
module HWBRAN and in the attendant related mod- 
ules [30] with the following substitutions: 



DGLAP-CS Pab 



IR-I DGLAP-CS PX"ji 



(3) 



while leaving the hard processes alone for the mo- 
ment. We have in progress [31]thc inclusion of YFS 
synthesized electroweak modules from Refs. [32]for 
HERWIG6.5, HERWIG-l-h [33] hard processes, as the 
CTEQ [34] and MRST(MSTW) [35] best (after 2007) 
parton densities do not include the precision elec- 
troweak higher order corrections that do enter in a 1% 
prccison tag biidgct for processes such as single heavy 
gauge boson production in the LHC environment [3] . 

For definiteness, let us illustrate the implementa- 
tion by an example [36, 37], which for pedagogical 
reasons we will take as a simple leading log shower 
component with a virtuality evolution variable, with 
the understanding that in HERWIG6.5 the shower de- 
velopment is angle ordered [36] so that the evolution 
variable is actually ~ E9 where 9 is the opening angle 
of the shower as defined in Ref. [36] for a parton ini- 
tial energy In this pedagogical example, which we 
take from Ref. [36], the probability that no branch- 
ing occurs above virtuality cutoff Qq is Aa{Q^, Qq) so 
that 

rfA„(i, Ql) = ^A(i, Ql) ^ J dz^P,a{z), (4) 



which implies 



Aa(<3^(3o) =exp 



(5) 

The attendant non-branching probability appearing in 
the evolution equation is 



^We thank M. Seymour and B. Webber for discussion on this 
point. 
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A((5^^) = ^f^^-T^f^, t = kl= virtuality of gluon a. 



Aa(t,Q?) 

The respective virtuality of parton a is then 
generated with 

AaiQ^t) = R, 



(6) 



(7) 



where ii is a random number uniformly distributed 
in [0, 1] . With 



asiQ) 



we get for example 



27r 



bo log (f) 



(8) 



/' 

Jo 



27r 



27r6oln(^) 



^^oln(f)- 
so that the subsequent integration over dt yields 

1 = 



(9) 




(10) 



Finally, introducing I into (5) yields 



K{Q\Qi) = exp 



In 



36o 



■In 



A2 



In 



(^) 



(11) 



If we now let Ao(Q^, i) = -R, then 



Mil 



3ET7 



R 



(12) 



which implies 



t = K' 



A2 



ij 2 



(13) 



Recall in HERWIG6.5 [11] we have 



(llric-lO), n/ = 5 



= -BETAF. 

3 



(14) 



where in the last line we used the notation in 
HERWIG6.5. The momentum available after a qq 
split in HERWIG6.5 [11] is given by 



QQBAR = QCDL,(^^^ 



(15) 



in complete agreement with (13) when we note the 
identifications 

t = QQBAR^, A = QCDL3, Q = QLST. 

The leading log exercise leads to the same algebraic 

relationship that HERWIG6.5 has between QQBAR 

and QLST but we stress that in HERWIG6.5 these 

quantities are the angle-ordered counterparts of the 

virtualities we used in our example, so that the 

shower is angle-ordered. 

When we repeat the above calculation for the 

IR-Improved kernels in (2), we have 



qG 

SO that 
dz- 



^2(l-^)TG + (l-z)2zT<5 

(16) 



/ 

Jo 



27r 



"""^^ ~ 6oln(f^)(7c + l)(70+2) 



(7G+3)- 

This leads to the following integral over dt 

4fyFs(7G)e^''° 



(17) 



dt 



Ql t 6oln(-^)(7G--^l)(7a+2)(7G+3) 
(io(7G + l)(7G+2)(7G+3) 



Ei 1 



8.369604402 \ 



Ql 



Here we have used 



asCc ( Tf_ 
3 



2 TT 

We finally get the IR-improved formula 

Aa{Q\t) = e^v[- {F {Q') - F {t))] , 
where 



Fm = 



bo (7G + 1) (7G + 2) (7G + 3) 
^. r 8.369604402 \ 



(18) 

(19) 

(20) 

(21) 
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and Ei is the exponential integral function. In 
Fig. 1 we show the difference between the two results 
for Aa(Q^,t). We see that they agree within a few % 
except for the softer values of t, as expected. We 
look forward to determining definitively whether the 
experimental data prefer one over the other. This 
detailed study will appear elsewhere [38] but we 
begin the discussion below with a view on recent 
FNAL data. Again, we note that the comparison in 




0.85 



0.8 



DGLAP-CS 
IR.Inip.DGLAP-CS 



4 6 
t (GeV^) 



10 



an oversight, as the module HWSFBR which 
controls the remainder of the space-like branching 
processes was also still in need of IR-improvement in 
versions 1.0 and 1.02"*. We have done the required 
IR-improvement of the latter module as well in 
version 1.031. While the effect in going from version 
1.0 to version 1.02 is small, that in going from 
version 1.02 to 1.031 is not in general and we 
recommend version 1.031 for the best precision. We 
now illustrate some of the results we have obtained 
in comparing ISR showers in IIERWIG6.5 and with 
those in HERWIRI1.0(31) at LHC and at FNAL 
energies, where some comparison with real data is 
also featured at the FNAL energy. Specifically, we 
compare the z-distributions, pr-distributions, etc., 
that result from the IR-improved and usual 
DGLAP-CS showers in what follows^. 
First, for the generic 2— >-2 hard processes at LHC 
energies (14 TeV) we get the comparison shown 
Figs. 2, 3 for the respective ISR z-distribution and 
distribution at the parton level. Here, there are 
no cuts placed on the MC data and we define z as 



parton 



/E^ 



where -Ebcam is the cms beam 



energy and E'parton is the respective parton energy in 
the cms system. The two quantities z and for 
partons are of course not directly observable but 
their distributions show the softening of the IR 
divergence as we expect. Turning next to the 



Energy Traction of parton shower constituents ^ QCD 2^2 hard parton scattering. 



Figure 1: Graph of Aa((^^^) for the DGLAP-CS and 
IR.Imp.DGLAP-CS kernels (11, 20). is a typical vir- 
tuahty closer to the squared scale of the hard sub-process 
- here we use — 25GeV^ for illustration. 

Fig. 1 is done here at the leading log virtuality level, 
but the sub-leading effects we have suppressed in 
discussing it will not change our general conclusions 
drawn therefrom. 

We have carried out the corresponding changes for 
all of the kernels in (2) in the HERWIG6.5 
environment, with its angle-ordered showers, 
resuhing in the new MC, HERWIRIl.O, in which the 
ISR parton showers have IR-improvement as given 
by the kernels in (3). In the original release of the 
program, v. 1.0, we stated that the time-like parton 
showers had been completely IR-improved in a way 
that suggested the space-like parton showers had not 
yet been IR-improved at all. Then, in release 1.02, 
we stated that the only part of the space-like parton 
showers without IR-improvement in v. 1.0 is that 
associated with HERWIG6.5's space-like module 
HWSGQQ for the space-like branching process 
G — 7> g?, a process which is not IR divergent and 
which is, in any case, a sub-dominant part of the 
shower. The module HWSQGG was IR-improved as 
weh in the release HERWIRI1.02. This was in fact 



DGLAP 
- IR.Imp.DCLAP 




0.001 



0.002 0.003 
Energy fraction 



0.004 



0.005 



Figure 2: The z-distribution(ISR parton energy fraction) 
shower comparison in HERWIG6.5. 

similar quantities for the tt^ production in the 



''We thank Profs. B. Webber and M. Seymour for discussion 
here. 

^Note that similar results for PYTHIA and MCONLO are in 
progress in general; for MCCSNLO we have some initial results 
already in particular cases — see the discussion below. 
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Figure 3: The P-^-distribution (ISR parton) shower com- 
parison in HERWIG6.5. 



Figure 5: The tt"^ Py-distribution shower comparison in 
HERWIG6.5. 



generic 2 — > 2 hard processes at LHC, we see in 
Figs. 4, 5 that spectra in the former are similar and 
spectra in the latter are again softer in the 
IR-improved case. These spectra of course would be 
subject to some "tuning" in a real experiment and 
we await with anticipation the outcome of such an 
effort in comparison to LHC data. 

Energy fraction of in HERWIG6.5 for QCD 2 -> 2 hard parton scattering. 



DGLAP 
- IR.Imp.DGLAP 




0.0001 



0.0002 0.0003 
Energy fraction 



0.0004 0.0005 



Z-rapidity distribution with cuts on the acceptance 
as 40GeV < Mz, Pt > 5GeV for Z ^ fifi - all 
lepton rapidities are included. For the energy 
fraction distribution and the pt distributions we 
again see softer spectra in the former and we see 
similar spectra in the latter in the IR-improved case. 
For the rapidity plot, we see the migration of some 
events to the higher values of \r]\, which is not 
inconsistent with a softer spectrum for the 
JR-improved case ^. We look forward to the 
'confrontation with experiment, where again we stress 
that in a real experiment, a certain amount of 
"tuning" with affect these results. The question will 
always be which set of distributions gives a better 
per degree of freedom. 

Finally, we turn the issue of the IR-cut-ofF in 
HERWIG6.5. In HERWIG6.5, there is are IR-cut-off 
parameters used to separate real and virtual effects 
and necessitated by the -I— function representation of 
the usual DGLAP-CS kernels. In HERWIRI, these 
parameters can be taken arbitrarily close to zero, as 
the IR-improved kernels are integrable [9, 10]. We 
now illustrate the difference in IR-cut-off response by 
comparing it for HERWIG6.5 and HERWIRI: we 
change the default values of the parameters in 
HERWIG6.5 by factors of .7 and 1.44 as shown in 
the Fig. 9. We see that the harder cut-off reduces 



Figure 4: The vr"*" energy fraction distribution shower com- 
parison in HERWIG6.5. 

We turn next to the luminosity process of single Z 
production at the LHC, where in Figs. 6,7,8 we show 
respectively the ISR parton energy fraction 
distribution, the Z-pT distribution, and the 



"One might wonder why we show the Z rapidity here as 
the soft gluons which we study only have an indirect affect 
on it via momentum conservation? But, this means that the 
rapidity predicted by the IR-improved showers should be close 
to that predicted by the un-improved showers and we show this 
cross-check is indeed fulfilled in our plots. 
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[Energy fraction distribution of parton shower for single Z production.[ 



The unit normalized differential cross section for Z production as a function of vector boson rapidity. 
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Figure 6: The 2:-distribution(ISR parton energy fraction) 
shower comparison in HERWIG6.5. 



Figure 8: The Z rapidity-distribution(ISR parton shower) 
comparison in HERWIG6.5. 
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Figure 7: The Z pT-distribution(ISR parton shower effect) 
comparison in HERWIG6.5. 



the phase space only significantly for the 
IR-improved kernels and that the softer cut-off has 
also a small effect on the usual kernels spectra 
whereas as expected the IR-improved kernels spectra 
move significantly toward softer values as a 
convergent integral would lead one to expect. This 
should lead to a better description of the soft 
radiation data at LHC. We await confrontation with 
experiment accordingly. 

We finish this initial comparison discussion by 
turning to the data from FNAL on the Z rapidity 
and pt spectra as reported in Refs. [39, 40]. We 
show these results, for 1.96TeV cms energy, in 



Fig. 10. We see that, in the case of the CDF rapidity 
data, HERWIRI1.0(31) and HERWIG6.5 both give a 
reasonable overall representation of the data but 
that HERWIRI1.0(31) is somewhat closer to the 
data for small values of Y . The two x^/d.o.f are 1.77 
and 1.54 for HERWIG6.5 and HERWIRI1.031 
respectively. The data errors in Fig. 10(a) do not 
include luminosity and PDF errors [39], so that they 
can only be used conditionally at this point. 
Including the NLO contributions to the hard process 
via MC@NLO/HERWIG6.510 and 
MC@NL0/HERWIRI1.G31[28]^ improves the 
agreement for both IIERWIG6.5 and for 
IIERWIRI1.031, where the x^/d.o.f are changed to 
1.40 and 1.42 respectively. 

For the DO pr data, we see that HERWIRI1.0(31) 
gives a better fit to the data compared to 
HERWIG6.5 for low pr, (for pt < 12.5GeV, the 
X^/d.o.f. are ^ 2.5 and 3.3 respectively if we add the 
statistical and systematic errors), showing that the 
IR-improvement makes a better representation of 
QCD in the soft regime for a given fixed order in 
perturbation theory. Including the results of 
MG@NLO [28] improves the xVd-O.f for the 
IIERWIRI1.031 in both the soft and hard regimes 
and it improves the HERWIG6.510 x^/d.o.f for px 
near 3.75 GeV where the distribution peaks. For 
PT < 7.5GeV the xVd-o.f for the 
MC@NLO/HERWIRI1.031 is 1.5 whereas that for 
MC@NLO/HERWIG6.510 is worse. 



^We thank S. Frixione for helpful discussions with this im- 
plementation. 
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Figure 9: IR-cut-ofT sensitivity in z-distributions of the ISR parton energy fraction: (a), DGLAP-CS (b), IR-I DGLAP-CS 
- for the single Z hard sub-process in HERWIG-6.5 environment. 
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Figure 10: Comparison with FNAL data: (a), CDF rapidity data on (Z/7*) production to e^e~ pairs, the circular 
dots are the data; (b), DO pr spectrum data on (Z/j*) production to e^e~ pairs, the circular dots are the data, 
the blue triangles are HERWIRI1.031, the green triangles are HERWIG6.510 - in both (a) and (b) the blue squares 
are MC@NLO/HERWIRI1.031, and the green squares are MC@NLO/HERWIG6.510. These are untuned theoretical 
results.. 



5. Conclusions 

Our new MC HERWIRI1.0(31) sets the stage for the 
further implementation of the attendant [4] new 
approach to precision QEDxQCD predictions for 
LHC physics by the introduction of the respective 
resummed residuals needed to systematically 
improve the precision tag to the 1% regime for such 
processes as single heavy gauge boson production, 
for example. Here, we already note that this new 
IR-improved MC, HERWIRI1.0(31), available at 



http://thep03.baylor.edu, is expected to allow for a 
better per degree of freedom in data analysis of 
high energy hadron-hadron scattering for soft 
radiative effects and we have given evidence that this 
is indeed the case. 
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